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In the random-interaction ensembles, three proportional correlations between quadrupole mo¬ 
ments of the first two N = 2~^ states robustly emerge, including Q(2^) = ±Q(2j) correlations 
consistently with realistic nuclear survey, and the Q(2^) = —1(3(2)^) correlation, which is only 
observed in the sd-boson space. These correlations can be microscopically characterized by the 
rotational SU(3) symmetry and quadrupole vibrational U(5) limit, respectively, according to the 
Elliott model and the sd-boson mean-field theory. The anharmonic vibration may be another phe¬ 
nomenological interpretation for the Q(2^) = —Q(2^) correlation, whose spectral evidence, however, 
is insufficient. 

PACS numbers: 21.10.Ky, 21.60.Cs, 21.60.Fw, 24.60.Lz 


I. INTERACTION 

Finite many-body systems (e.g., nuclei, small metal¬ 
lic grains, metallic clusters) robustly maintain similar 
regularities, despite their different binding interactions. 
For example, they all present the odd-even staggering on 
their binding energies, which are, however, attributed to 
various mechanisms [iH^- Particularly in nuclear sys¬ 
tems, the nucleon-nucleon interactions numerically ex¬ 
hibit a “random” pattern with no trace of symmetry 
groups, whereas nuclear spectra follow some robust dy¬ 
namical features: the nuclear spectral fluctuation is uni¬ 
versally observed [1-Q; low-lying spectra of even-even 
nuclei are orderly and systematically characterized by se¬ 
niority, vibrational and rotational structures 00 , be¬ 
yond = O'*' ground states without exception. 

To demonstrate the insensitivity of these robust regu¬ 
larities to the interaction details, and to reveal its under¬ 
lying origin, random interactions are employed to simu¬ 
late (or even introduce) the variety and chaos into a finite 
many-body system. Thus, the predominant behaviors in 
a random-interaction ensemble correspond to dynamical 
features in a realistic system. Many efforts have been 
devoted along this direction fTTI - [l5l| . For instance, sim¬ 
ilarly to realistic even-even nuclei, the predominance of 
the 7 = 0 ground states [il,[i3 and collective band struc¬ 
tures [H, [I^ have been observed in random-interaction 
ensembles. However, there are only few attempts to 
study the robustness of nuclear quadrupole collectivity 
against the random interaction. This is partly because 
a random-interaction ensemble potentially gives weaker 
E2 transitions than a shell-model calculation with “re¬ 
alistic” interactions [1^. Even so, some robust corre¬ 
lations about the E2 collectivity can be expected. For 
example, the Alaga ratio between the quadrupole mo¬ 
ment (Q) of the 2^ state and B(E2, 2f —s- 0)*') high- 
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FIG. 1. (Color online) (2+||E 2||2+ ) and (2+||E2||2+) matrix 
elements from Table I of ref. (be., [Q{2^), Q(2j)] plots 
scaled by ^167r/5(2220|22)). The 6 parameterization defined 
in Eq. © is illustrated in the red sector. The Q(2^) = 
—Q{2f) correlation is obvious along the 6 = —45° direction 
(the black dash diagonal line). 


lights both near-spherical shape and well deformed rotor 
in random-interaction ensembles [H, HH; ratios of E2 
transition rates between yrast 0)^, 2+ and 4^ states are 
also correlated to the ratio of 2'^ and 4^ excitation ener¬ 
gies 

This work further studies the robust correlation be¬ 
tween Q values of the hrst two 2^ states, inspired by a 
recent experimental survey [2^. As shown in Fig. [U this 
survey demonstrated a global Q(2^) = —Q{2^) corre¬ 
lation across a wide range of masses, deformations, and 
2^ energies. We will make use of random-interaction en- 
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sembles to provide an interacting-particle vision to this 
correlation, and search for other underlying Q correla¬ 
tions. The statistic analysis based on the Elliott SU(3) 
model [2^ and the mean-field Hartree-Bose theory 
is applied. 


II. CALCULATION FRAMEWORK 


In our random-interaction calculations, the single- 
particle-energy degree of freedom is switched off to avoid 
the interference from the shell-structure detail. The two- 
body interaction matrix element, on the other hand, 
is denoted by usual, where ji, j2, js and 

j4 represent the angular momenta of single-particle or¬ 
bits (half integer for fermions and integer for bosons), 
and the superscript J labels the total angular momen¬ 
tum of the two-body configurations involved the inter¬ 
action element. In our calculations, is random¬ 
ized independently and Gaussianly with (/r = 0, = 

1 which insures the invariance of our ran¬ 

dom two-body interactions under arbitrary orthogonal 
transformations [2^. All the possibilities of random in¬ 
teractions and their outputs via microscope-calculations 
construct the two-body random ensemble (TBRE) [23- 
[2^. Obviously, in the TBRE, diagonal interaction ele¬ 
ments potentially have larger magnitudes. 

For the shell-model TBRE in this work, four model 
spaces with either four or six valence protons in either sd 
or pf shell are considered, correspondingly to four nuclei: 
24gi, 26s^ 44(^J. 46pg TBRE, sd-boSOn 

spaces are constructed for nuclei with valence boson num¬ 
bers Nb =12, 13, 14 and 15, where s and d represents 
I = Oh and I = 2h bosons, respectively. It is noteworthy 
that a single calculation with random interactions does 
not match, and does not intend to match, to a realistic 
nucleus. It only presents a pseudo nucleus in the compu¬ 
tational laboratory. Thus, in this article, model spaces 
described above are named as corresponding pseudo nu¬ 
clei for convenience. For example, the model space with 
four protons in the sd shell corresponds to pseudo ^^Si. 
Statistic properties of many random-interaction calcula¬ 
tions for pseudo nuclei can be related to the robustness of 
dynamic features in realistic nuclear systems. To insure 
the statistic validity of our conclusions, 1 000 000 sets of 
random interactions are generated for each pseudo nu¬ 
cleus, and inputted into the shell-model or IBMl calcu¬ 
lations. If one calculation produces a / = 0 ground state, 
Q matrix elements of 2'^ and 2^ states will be calculated 
and recorded for the following statistic analysis. 


III. Q CORRELATIONS IN SHELL MODEL 

In the Shell Model, the Q matrix element of one 2+ 
state, |2+), is defined conventionally as 

Q(2+) = {2+||Q||2+), 

Q = {j\\r^Y^\\j'){al X 

where a| and dji are single-particle creation and time- 
reversal operators at orbits j and j', respectively. A pro¬ 
portional Q correlation between the first two 2f states 
is normally characterized by the ratio of (5(2j)/Q(2)''). 
Geometrically, such correlation also corresponds to a 
straight line with the polar angle, 

0 = arctan|^^l, (2) 

across the origin in the [Q{2f), <5(2^)] plane. For ex¬ 
ample, the experimental Q{2'^) = —(5(2^) correlation 
suggested by Ref. can be illustrated by a diago¬ 
nal 9 = arctan(—1) = —45° line as expected in Fig. [T] 
We also visualize the polar-angle scheme of the [(5(2)*"), 
(5(2J)] plane in Fig. [J 

In this work, we prefer the statistic analysis based 
on the polar angle 9 over the Q{22)/Q{2f) ratio be¬ 
cause of two reasons. Firstly, the distribution of the 
Q{ 22 )/Q{ 2 '^) ratio spreads widely, so that the statis¬ 
tic detail about ( 5 ( 2 ^) = —Q{2^) correlation may be 
concealed. In particular, there robustly exists 8% prob¬ 
ability of \Q{ 22 )/Q[2^)\ > 10 due to the predominance 
of weak quadruple collectivity, i.e. small |(5(2('')|, in the 
shell-model TBRE [2^. However, we intend as compre¬ 
hensively as possible to present the statistic detail about 
the experimental (5(2^) = —Q{2f) = —1 correlation. 
The wide statistic range of the (5(2j)/(5(2('’) ratio may 
fails this intention. By converting the (5(2j)/(5(2)'") ra¬ 
tio to the 9 value, the statistic range is limited between 
(—90°, 90°), and a clearer vision around the (5(2^) = 
—Q{2^) correlation can be obtained around 9 = —45°. 
Secondly, the 9 parameterization intuitively provides a 
reasonable geometric standard of symmetric sampling. 
Taking the (5(2 J) = —Q{2^) correlation for example, 
there is actually no (pseudo) nucleus following exact 
(5(2^) = —Q{2^) relation in experiments or our TBRE, 
and yet we can take 9 G (—50°,—40°) as the sampling 
range to represent this correlation. One sees this sam¬ 
pling range indeed covers a symmetric area related to 
the (5(2^) = —Q{2f) correlation in the [Q{2f), (5(2^)] 
plane. With the Q{22)/Q(2f) statistic, the determina¬ 
tion of symmetric sampling range for one specifically Q 
correlation can be controversial or simply another repre¬ 
sentation of the 9 parameterization. Therefore, all the 
statistic, analyses, and discussions in this work are based 
on the 0 value. 

In Fig. [21 we present 9 distributions of four pseudo 
nuclei in the shell-model TBRE compared with the ex¬ 
perimental distribution from Ref. |23l ]. The experimen¬ 
tal Q{2^) = —(5(2^) correlation is represented by the 
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FIG. 2. (Color online) 6 distributions from the experimen¬ 
tal survey (Exp) [23 ] and the shell-model TBRE. 9 = ±45° 
peaks are highlighted correspondingly to <5(2J) = ±Q{2^) 
correlations, respectively. Error bars correspond to statistic 
error. 


main peak around 6 = —45°, which is also reproduced 
by the TBRE. Furthermore, several weak peaks around 
6 = 45° are also observed in both experimental data 
and random-interaction systems, corresponding to the 
( 5(22 ) = Q{‘2t) correlation. 

As proposed by Ref. [l^, nuclear rotor models can 
give the 0 = —45° correlation, even although such cor¬ 
relation experimentally occurs in both rotational or non- 
rotational nuclei. Therefore, we will further examine 
whether 0 = ±45° correlations is the symbol of the un¬ 
derlying rotational collectivity in TBRE. Firstly, we ver¬ 
ify whether 0 = ±45° correlations accompany rotational 
spectra in the TBRE. Secondly, we search statistic sig¬ 
nature of the random-interaction elements that provides 
the 9 = ±45° correlations, and trace such signature back 
to the microscopical Hamiltonian of nuclear rotor model, 
namely the Elliott SU(3) Hamiltonian. 

Following previous random-interaction studies [H, [l^, 
HI) 113) potential rotational spectra with 9 = ±45° 
correlations can be characterized by the energy ratio 
i ?42 = E^+/E^+ ~ 10/3, where E^+ and E^+ corre¬ 
spond to the excitation energy of yrast 2+ and 4+ states, 
respectively. Thus, we plot the R 42 distributions with 
9 G (—50°, —40°) and 9 G (40°, 50°), respectively, in Fig. 
[31 and compare them with that in the whole TBRE. Ex¬ 
cept for R distributions in both 9 G (—50°,—40°) 
and 9 G (40°, 50°) regions are identical to those in the 
whole TBRE within statistic error. For the R 42 dis¬ 
tribution has an observable enhancement at R 42 = 1 with 
9 G (40°, 50°). Namely, the 9 = 45° correlation seems to 
partially originate from the seniority-like level scheme in 
space. This observation explains why the 9 = 45° 
peak for is stronger shown in Fig. |31 given the domi¬ 
nance of pairing-like behaviors in the TBRE [ilElIIl. 
Nevertheless, there is no special favor on rotational spec¬ 
tra from 9 = ±45° correlations in the shell-model TBRE, 
consistently with the survey on the realistic nuclear sys- 
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FIG. 3. (Golor online) R 42 distributions around 6 = ±45° 
correlations (red circles and blue triangles, respectively) com¬ 
pared with those in the whole shell-model TBRE (black 
squares). Error bars correspond to statistic error. 


tern [23(. 

In Ref. 2l|, the interaction signature of prolate and 
oblate shapes is represented by the average values of in¬ 
teraction elements (denoted by In this work, 

we also adopt V/fto probe the the interaction sig- 
nature oi 9 = ±45° correlations. In detail, we collect 
all the interaction elements within 9 G (—50°,—40°) 
and 9 G (40°, 50°), normalize them by the factor of 

E Vnnj.u ’ then calculate all the val- 

Jjihhh 


ues for both 9 G (—50°, —40°) and 9 G (40°, 50°) regions, 
respectively. Because signs of interaction elements can be 
changed by different phase conventions, we only discuss 
magnitudes of (denoted by to avoid 

the potential ambiguity from phase conventions. To sim¬ 


plify the following discussion, each I labeled by 

the index, 2ji x 10000±2j2X 1000±2/3X 100±2j4X 10± J. 
For example, the pairing force between S 4/2 or pi /2 nu¬ 
cleons, Eiiii, corresponds to index “11110”. We list 
_ 2 2 2 2 


IE/,- , , I values of both 9 G (-50°,-40°) and 9 G 

J L J 2 J O J ^ 

(40°, 50°) region in an increasing order of their indices 
in Tables III and ITll 


To comprehensively compare II// 2 J 3 J 4 I values between 
9 G (—50°,—40°) and 9 G (40°, 50°) regions, we plot 
them against their order numbers (see Table U and [TTI 
in Fig. m Most of E/, . . are close to zero following 
the ensemble distribution. However, there are several 
relatively large |E/ | values, which presents obvious 

peaks in Fig. |4l Peak positions for 9 G (—50°, —40°) are 
roughly consistent with those for 9 G (40°, 50°), which 
hints that 9 = ±45° correlations may share the same 
interaction signature. 
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FIG. 4. (Color online) l(0ii2)‘^iC'su(3)l(i3j4)'^)l values (see text for definitions) against order numbers from 

Tables U and im indices are highlighted for obvious peaks for |((jij 2 )'^|C'su( 3 )l(j 3 j 4 )‘^)| values. 


The interaction signature of 0 = ±45° correlations can 
be related to the Elliott Hamiltonian. Such Hamiltonian 
is dominated by the SU(3) Casimir operator as defined 
by 

^su( 3 ) = -^Q ■ Q + ■ L, (3) 

where Q and L are quadrupole-moment and orbital- 
angular-momentum operators. We calculate matrix ele¬ 
ments of ((jij 2 )'^|(?su( 3 )l(j 3 j 4 )'^)! and still focus on their 
magnitudes (denoted by |(jij 2 |C'su( 3 )|j 3 j 4 )‘^|), similarly 

to the treatment for Kjij 2 |C'su( 3 )|j 3 j 4 )‘^| is also 

labeled by index, 2ji x 10000 ± 2 j 2 x 1000 ± 2j^ x 100 ± 
2 j 4 X 10 ± J, and thus comparable with 
shown in Table mini and Fig. |4l In Fig. |41 relatively 
large Kjij 2 |f?su( 3 ) Ij 3 j 4 )‘^l also presents several obvious 
peaks, which have similar pattern to peaks for 

both e G (-50°,-40°) and 9 G (40°, 50°) regions. This 
observation implies the relation between the SU(3) sym¬ 
metry and 0 = ±45° correlations. 

We also highlight indices for | (jij 2 1Csufs) \hk) ^ I peaks 
in Fig. 01 according to which, the SU(3) Casimir operator 
always has large magnitudes for diagonal matrix elements 
with jij 2 = 33 j A- On the other hand, large \ for 

9 = ±45° correlations also occurs for diagonal ji j 2 = jsjA 
in Table m and HH As described in Sec. m larger magni¬ 
tudes of diagonal elements is required by the invariance 


of TBRE under orthogonal transformation of two-body 
configuration. Therefore, the shell-model TBRE intrinsi¬ 
cally maintains part of the SU(3) properties to restore the 
9 = ±45° correlations, even though it spectrally presents 
no trace of the SU(3) symmetry as illustrated in Fig. 01 

After clarifying the relation between 9 = ±45° cor¬ 
relations and the SU(3) symmetry, we microscopically 
describe how these two Q correlations emerge in a major 
shell, i.e. sd or pf shell here. In the Elliott model, any 
2+ state within a major shell is labeled by the SU(3) rep¬ 
resentation (A, p), the quantum number of the intrinsic 
state (AT), and orbital angular momentum L = 2 [^ . 
The 2+ state is normally near the bottom of a K band, 
and thus its Q value can be approximately given by [s^ 

p \ 

Q{2+) = -{K^-2). (4) 

The K number is limited to 0, I and 2. Thus, the 
Q{2.f) = —(5(2^), i.e. 9 = —45°, correlation is produced 
by two 2+ states with the same A number and AT = 0, 2 
respectively, which agrees with the rotor-model conjec¬ 
ture [1^. On the other hand, the Q{2t) = i-®- 

9 = 45°, correlation is from two 2+ states with the same 
A and K values. 

According to above SU(3) description, one can expect 
two 2+ states with the 9 = —45° correlation from the 
same (A, p) representation. On the contrary, a single 
(A, p) representation can not produced two 2“'' states 
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TABLE I. values around 6 ±45° correlations and 

IOij2|C'su(3)|j3i4)'^| elements [“SU(3)” column] in the sd 
shell. The “index” column presents the integer, 2ji x 10000 ± 
2 j 2 X 1000 ±2j3 X 100 ±2^4 X 10 ±J, to identify two-body in¬ 
teraction elements. All the data is organized in an increasing 
order of the index column. 


Order 

Index 

1 ^71 .;2.73.;4 1 

SU(3) 

e G (-50°,- 

-40°) e G (40°, 50°) 

1 

11110 

0.031 

0.028 

20.0 

2 

11330 

0.009 

0.006 

5.7 

3 

11550 

0.017 

0.006 

6.9 

4 

13131 

0.047 

0.007 

7.0 

5 

13132 

0.022 

0.053 

10.2 

6 

13152 

0.007 

0.028 

3.9 

7 

13332 

0.035 

0.025 

2.5 

8 

13351 

0.006 

0.007 

0.0 

9 

13352 

0.002 

0.043 

2.3 

10 

13552 

0.006 

0.024 

3.3 

11 

15152 

0.045 

0.087 

11.8 

12 

15153 

0.095 

0.172 

7.0 

13 

15332 

0.024 

0.005 

3.1 

14 

15352 

0.016 

0.014 

2.8 

15 

15353 

0.004 

0.060 

0.0 

16 

15552 

0.064 

0.002 

4.0 

17 

33330 

0.118 

0.039 

15.0 

18 

33332 

0.136 

0.107 

0.7 

19 

33352 

0.025 

0.021 

4.1 

20 

33550 

0.014 

0.014 

2.4 

21 

33552 

0.016 

0.018 

0.4 

22 

35351 

0.018 

0.076 

13.0 

23 

35352 

0.015 

0.104 

7.2 

24 

35353 

0.180 

0.009 

2.0 

25 

35354 

0.028 

0.167 

2.0 

26 

35552 

0.015 

0.028 

4.8 

27 

35554 

0.046 

0.021 

0.0 

28 

55550 

0.151 

0.086 

16.0 

29 

55552 

0.046 

0.069 

8.1 

30 

55554 

0.160 

0.037 

2.0 

with the 

same . 

K number, so that the 9 = 45° 

corre- 


lation always requires the cooperation of two different 
(A, /r) representations. Empirically, the former case has 
a relatively larger probability to emerge in the low-lying 
region, which explains why the 9 = —45° peak intensity 
is always larger than the 9 = 45° one in Fig. [2j 

Independently of the rotor interpretation, the anhar- 
monic vibration (AHV) with quadrupole degrees of free¬ 
dom [s^ can also provide the 9 = —45° correlation. In 
the AHV interpretation, the first two 2+ states are con¬ 
structed with a significant mixing of one- and two-phonon 
configurations as 

\2t)=a,\b^) + a2\{b^f), 

|2+) = -a2|6t)±ai|(6t)2)^ 


where b^ is the creation operator of a phonon; oi and 02 
are amplitudes of phonon configurations. In this phonon 
space, the quadrupole operator Q is a polynomial of op¬ 
erator fel ± 5 (^ . where b is the phonon time-reversal 
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FIG. 5. (Color online) R 22 distributions around 6 — —45° 
(red circles) compared with those in the whole shell-model 
TBRE (black squares). Error bars correspond to statistic 
error. 


operator. The first order of such polynomial dominates 
the Q matrix element. However, it also vanishes with re¬ 
spect to conhguration with definite numbers of phonons. 
In particular, 

(6||Q||6^) oc (6||5t ±5||6^) = 0, 

((6)2||Q||(5t)2)oc((5)2||6t±6||(5t)2) = o, 


Thus, 


(2+||Q||2+)=2aia2(6||Q||(6t)2), 

(2+||Q||2+)=-2aia2(5||Q||(6t)2), 


and the 9 = —45° relation is obtained. 

We can spectrally examine this AHV interpretation 
for the 9 = —45° correlation in the shell-model TBRE. 
Because AHV 2+ states correspond to the mixing of one- 
and two-phonon conhgurations as defined in Eq. the 
excitation energy of the first 2+state, E{2f), is smaller 
than the one-phonon excitation energy, feu; while i?(2^) 
is larger than 2Iku, according to the perturbation theory. 
Thus, the energy ratio of R 22 = A(2j)/E(2]'') of the 
AHV is always larger than 2. In other words, if the AHV 
contributes to the 9 = —45° correlation in the TBRE, 
the distribution of R 22 with 9 € (—50°,—40°) should 
have an obvious enhancement for R 22 > 2. In Fig. [5l we 
compare R 22 distributions in the 9 G (—50°, —40°) range 
and those in the whole shell-model TBRE. There is no 
obvious difference between these R 22 distributions. Thus, 
we don’t see the spectral sign of the AHV contribution 
to the 9 = —45° correlation. 
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TABLE II. The same as Table U except for the pf shell. 


Order 

Index 

\V^ ■ - ■ 1 

_1 3132333A ' _ 

SU(3) 

Order 

Index 

\VJ. . . 1 

_1 J1J2J3J4 1_ 

SU(3) 

0 e (-50b-40“; 

) 9 G (40“,50“) 

0 G (-50“,-40“; 

) 6 G (40“,50“) 

1 

11110 

0.030 

0.011 

32.0 

48 

35352 

0.059 

0.110 

18.9 

2 

11330 

0.011 

0.008 

12.2 

49 

35353 

0.044 

0.030 

9.1 

3 

11550 

0.010 

0.006 

9.7 

50 

35354 

0.044 

0.011 

12.3 

4 

11770 

0.006 

0.007 

0.0 

51 

35372 

0.002 

0.001 

2.3 

5 

13131 

0.018 

0.033 

23.4 

52 

35373 

0.004 

0.005 

3.8 

6 

13132 

0.039 

0.004 

27.9 

53 

35374 

0.006 

0.007 

3.8 

7 

13152 

0.000 

0.000 

6.6 

54 

35552 

0.003 

0.011 

3.7 

8 

13332 

0.009 

0.007 

3.2 

55 

35554 

0.010 

0.000 

3.3 

9 

13351 

0.000 

0.000 

0.0 

56 

35571 

0.006 

0.000 

0.0 

10 

13352 

0.002 

0.000 

3.5 

57 

35572 

0.002 

0.002 

4.1 

11 

13372 

0.007 

0.004 

8.6 

58 

35573 

0.010 

0.004 

6.2 

12 

13552 

0.003 

0.005 

4.1 

59 

35574 

0.005 

0.002 

5.1 

13 

13571 

0.008 

0.000 

9.0 

60 

35772 

0.004 

0.001 

3.3 

14 

13572 

0.007 

0.005 

5.4 

61 

35774 

0.004 

0.005 

2.6 

15 

13772 

0.004 

0.002 

0.0 

62 

37372 

0.035 

0.032 

26.8 

16 

15152 

0.006 

0.006 

20.1 

63 

37373 

0.033 

0.041 

15.3 

17 

15153 

0.077 

0.079 

6.5 

64 

37374 

0.067 

0.058 

12.9 

18 

15173 

0.002 

0.003 

1.0 

65 

37375 

0.060 

0.042 

6.0 

19 

15332 

0.001 

0.000 

4.7 

66 

37552 

0.001 

0.004 

0.9 

20 

15352 

0.005 

0.005 

7.9 

67 

37554 

0.000 

0.008 

1.0 

21 

15353 

0.001 

0.009 

1.3 

68 

37572 

0.004 

0.002 

2.4 

22 

15372 

0.000 

0.001 

2.8 

69 

37573 

0.001 

0.002 

3.3 

23 

15373 

0.000 

0.009 

1.9 

70 

37574 

0.007 

0.000 

1.5 

24 

15552 

0.025 

0.000 

6.1 

71 

37575 

0.005 

0.004 

0.0 

25 

15572 

0.003 

0.000 

2.4 

72 

37772 

0.010 

0.003 

3.4 

26 

15573 

0.008 

0.000 

2.2 

73 

37774 

0.026 

0.000 

8.6 

27 

15772 

0.001 

0.000 

0.0 

74 

55550 

0.071 

0.073 

26.4 

28 

17173 

0.069 

0.097 

18.9 

75 

55552 

0.041 

0.040 

12.2 

29 

17174 

0.045 

0.074 

12.0 

76 

55554 

0.042 

0.019 

7.1 

30 

17353 

0.000 

0.002 

1.1 

77 

55572 

0.007 

0.004 

5.8 

31 

17354 

0.011 

0.005 

5.9 

78 

55574 

0.003 

0.000 

3.6 

32 

17373 

0.014 

0.003 

8.9 

79 

55770 

0.016 

0.013 

2.1 

33 

17374 

0.015 

0.002 

2.0 

80 

55772 

0.009 

0.005 

1.2 

34 

17554 

0.002 

0.002 

3.4 

81 

55774 

0.002 

0.010 

0.2 

35 

17573 

0.010 

0.002 

4.8 

82 

57571 

0.006 

0.026 

24.6 

36 

17574 

0.013 

0.000 

6.3 

83 

57572 

0.067 

0.080 

18.0 

37 

17774 

0.000 

0.002 

0.0 

84 

57573 

0.053 

0.015 

7.2 

38 

33330 

0.070 

0.045 

40.6 

85 

57574 

0.028 

0.004 

0.2 

39 

33332 

0.101 

0.126 

25.6 

86 

57575 

0.095 

0.045 

9.0 

40 

33352 

0.005 

0.005 

2.5 

87 

57576 

0.020 

0.036 

9.0 

41 

33372 

0.007 

0.001 

6.1 

88 

57772 

0.010 

0.006 

5.4 

42 

33550 

0.011 

0.021 

2.0 

89 

57774 

0.010 

0.000 

5.1 

43 

33552 

0.010 

0.010 

0.2 

90 

57776 

0.013 

0.002 

0.0 

44 

33572 

0.004 

0.001 

4.4 

91 

77770 

0.082 

0.071 

27.0 

45 

33770 

0.009 

0.019 

10.2 

92 

77772 

0.025 

0.011 

18.6 

46 

33772 

0.005 

0.004 

6.7 

93 

77774 

0.092 

0.080 

3.1 

47 

35351 

0.004 

0.027 

27.0 

94 

77776 

0.123 

0.118 

9.0 


IV. g CORRELATIONS IN IBMl 


9 coordinates as 


In IBMl, the Q operator is a linear combination of two 
independent rank-two operators as: 


Q — + xQ^ 1 ( 8 ) 


9^ = arctan 


= arctan 


(2+||gi||2+) ] 

(2+||Qi||2+)J 

(2^||Q^||2+)J 


(9) 


where s + sd, Q'^ = [d^d]'^, and x is a free parame¬ 

ter. Correspondingly, we need to define two independent 


A robust correlation with the polar angle 9 should be 
insensitive to the x value, which requires = 02 = d. 
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Obviously, such correlation corresponds to a peak at (0, 
9) point in the two-dimensional {9^, 9^) distribution of 
the IBMl TBRE. 

Fig. [S] represents {9^, 9^) distributions of the IBMl 
TBRE with Nb = 12, 13, 14 and 15. These distribu¬ 


tions follow similar pattern with three sharp peaks along 
the 9^ = 9^ diagonal line, corresponding to three propor¬ 
tional Q correlations. We fit {9^, 0^) distributions to a 
two-dimensional function, f{9^, 9^), with three Gaussian 
peaks as 


3 

/( 6 »\ 6 » 2 ) ^ fg + '^A,exp 
2=1 


[(01 


0G) coswi -I- (0^ - 



where /o is the background; all the other fitting variables 
are parameters of Gaussian peaks. These three Gaussian 
peaks are labeled by indices z = 1, 2 and 3. For the 
zth peak, uJi defines its orientation in the (0i,0^) plane, 
(0G,0^j) is the peak position, Ai is the amplitude, and 
are widths along and perpendicularly to uji 
direction. Thus, the best-fit intensity of the zth peak can 
be calculated as 27r^iW|| 

In Table [ml we list the best-fit peak positions and 
intensities for the three sharp peaks in Fig. jS] The i = 
1 and i = 2 peaks are very close to (±45°, ±45°), i.e. 
“(5(2^) = ±Q{2^)” correlations, as labeled in Fig. |l]and 
Table iroi The i = 3 peak is located around (—21°, —22°), 
and thus gives Q{22)/Q{2f) ~ —3/7, the typical IBMl 
Q ratio at the U(5) limit regardless of the boson number. 
Therefore, we believe the z = 3 peak may correspond to 
the vibrational U(5) collectivity, and denote it as “U(5)” 
in following analysis. 

To identify or confirm the collective patterns corre¬ 
sponding to the three sharp peaks in Fig. | 6 l we firstly 
investigate their i ?42 distribution, i.e. the predominance 
of low-lying collective excitations, similarly to our i ?42 
analysis for the shell-model TBRE with Fig. [3] secondly, 
we adopt the sd-boson mean-field theory to observe dom¬ 
inant nuclear sharps of these peaks. 

For the analysis of Ra 2 distributions, we firstly collect 
all the random interactions, which produce (0^, 0^) points 
within 3° from peaks in Fig. [51 Secondly, all i ?42 values 
from these interactions are calculated. Thirdly, i ?42 dis¬ 
tributions of these peaks are calculated and presented 
in Fig. [71 respectively. <5(2^) = ±Q{2^) peaks always 
have large probabilities at rotational limit i ?42 = 3.3, 
which agrees with the rotor-model description. On the 
other hand, i ?42 distributions of U(5) peaks are domi¬ 
nated by i ?42 = 2, corresponding to a typical U(5) vibra¬ 
tional spectrum, which supports our U(5) assignment for 
this peak. 

Our analysis with the sd-boson mean-field theory 
starts with the sd-boson coherent state for the ground 
band as 


|g) = (st ±tanaodj)^'’|0). (11) 

Similarly to Ref. [s^, the nuclear shape, i.e. the opti¬ 
mized ao value, is determined by minimizing the Hamil- 


0/ J sinwi]2 [-(0i - 0/J siiioji ± (0^ - 0^ J cos 




( 10 ) 


tonian expectation value of this coherent state as 
Egla) =ai sin"^ a + a 2 sin^ a cos a 

^ . 2 2 ^ 4 ( 12 ) 

± 03 sm a cos a ± oq cos a, 

where Eg{ao) reaches the minimum of this equation; oq, 
ai, 02 and 03 are linear combination of sd-boson two- 
body interaction matrix elements as formulated in Ref. 
[ 3 ^. We calculate ao values for all the interactions with 
spin-Oli ground states in the TBRE, and perform fre¬ 
quency counting for calculated ao values. Thus, the 
ensemble-normalized a distribution for the zth peaks is 
given by 


P{a) = N{a,9E,9lA/Af{a) 


(13) 


where N(a,9E,9P) is the counting number with ao S 
(a-2.5°,a±2.5°) and ^(0i - 0±)2 + (02 - 02 < 30 ^ 

and A/"(a) is that with ao G (a — 2.5°, a ± 2.5°) in the 
whole IBMl TBRE. 

Fig. | 8 ] presents calculated P(a)s. The U(5) peak only 
has a significant probability at a = 0, corresponding 
to the s-boson condensation. Thus, 2+ states for the 
U(5) peak are constructed by replacing s bosons with d 
bosons in the s-boson condensation, which agrees with 
the quadrupole vibration described by the U(5) limit. 
This further confirms our U(5) characterization of this 
peak. On the other hand, the Q(2^) = ±Q(2/) peaks 
both have large probabilities for 0 < |a| < 90°, corre¬ 
sponding to the axially symmetric rotor at the SU(3) 
limit. Considering that the <5(2^) = ±Q{2^) peaks also 
favor SU(3) rotational spectra with R = 3.3 in Fig. 13 
we conclude that both Q( 2 j) = ±Q( 2 /) correlations in 
IBMl are strongly related to the SU(3) limit. 

Conversely, we also derive Q(2j) = ±Q(2/) correla¬ 
tions from the SU(3) limit of the IBMl. At the SU(3) 
limit, the 2f state is from the ground band with (A = 
2Nb, /z = 0) and K = 0\ yet, the 2^ state belongs to the 
(A = 2Nb — 4, ji = 2) representation, which gene rates /3 
and 7 bands with K = 0 and 2, respectively [2^. Thus, 
the 2 j is from either /3 or 7 band, which leads to two 
phase-different Q correlations 


Q^{2+) 4Nb-i Q^{2+) 4Nb-3 

Q9(2+) “ 4iVh±3’ Qa{2+) ~ 4iVh±3' 


( 14 ) 














FIG. 6. (Color online) Two-dimensional (6^, 6^) distributions of the IBMl TBRE. Three sharp peaks are characterized with 
“(3(2^) = ±(3(2)*")” and “U(5)” corrlations. 


TABLE III. Best-fit peak positions {9l i, 6^ ^) and intensities of three sharp peaks in Fig. [6]with the two-dimensional three-peak 
Gaussian function defined in Eq. (nsi). 



Qin) = - 

Q{n) 

Q{‘2t) = 

Q(2r) 

U(5) 


Nb 


Intensity 


Intensity 

9c, 3 ^c,3 

Intensity 


(deg) (deg) 

(xlO^ counts) 

(deg) (deg) 

(xlO^ counts) 

(deg) (deg) 

(xlO^ counts) 

12 

-42.05(1) -36.27(2) 

361(6) 

40.19(2) 43.44(1) 

272(5) 

-21.27(4) -22.17(1) 

565(9) 

13 

-42.29(1) -36.92(2) 

476(6) 

40.60(1) 43.62(1) 

347(5) 

-21.12(3) -22.21(1) 

742(9) 

14 

-42.52(1) -37.45(1) 

468(6) 

40.91(1) 43.73(1) 

347(5) 

-21.20(2) -22.23(1) 

733(8) 

15 

-42.66(1) -37.94(1) 

393(7) 

41.20(1) 43.83(1) 

334(5) 

-21.32(2) -22.26(1) 

683(8) 


For Nb oo, (3'^(2+) = Qs(2+) and Q^{2+) = -Qs(2+) 
are achieved, corresponding to <3(2^) = Q{2^) and 
( 3(22 ) = —Q{2f) correlations, respectively. 

In the shell-model TBRE, the ( 3(2 J) = —Q{2f) corre¬ 
lation has a larger probability than ( 3(2 J) = Q{2f) one 
(see Fig. [2]). Yet, in the IBMI TBRE, these two cor¬ 
relations have roughly equal peak intensities, i.e. prob¬ 
abilities, as shown in Fig. |6] and Table IIIIl This is a 
major difference between behaviors of the Q correlations 
in shell-model and IBMI TBREs. This difference can 
be explained according to the (A, fi) assignment of the 
SU(3) scheme. In the Shell Model, i.e. the Elliott model, 
the ( 3(2 J) = —Q(2f) correlation normally emerges with 


two 2+ states from a single (A, /r) representation, which 
empirically provides a larger probability. However, in 
the sd-boson space, both ( 3(22 ) = ±<3(2)^) correlations 
require 2^ and 2^ states from two different (A, /r) repre¬ 
sentations, and thus have similar probabilities. 

As shown in Table IIIIl J and J of (3(2^) = 
±Q{2^) correlations are systematically smaller than 45°. 
This observation can be explained with Eq. (fHl) . For 
large but finite Nb, the magnitude of (3(2^) is always 
smaller that that of Q(2f), which drives the \0N\ and 
\0lj\ value smaller than 45°. Therefore, we attribute the 
systematical derivation of (3(2^) — ±<3(2)'’) peak posi¬ 
tions from the exact SU(3) prediction to the finite-boson- 
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FIG. 7. (Color online) R 42 distributions around three peaks 
in Fig. [6] 



V. SUMMARY 


To summarize, we observe three proportional correla¬ 
tions between Q values of the first two I'" = 2+ states 
in the TBRE. = ±<5(2^) correlations robustly 

and universally exists in both shell-model and sd spaces, 
consistently with experiments. In the IBMl TBRE, the 
Q(2j) = correlation is also reported. By us¬ 

ing the Elliot model and the sd-boson mean-field the¬ 
ory, we can microscopically assign Q{2i) = ±(5(2J) 
correlations to the rotational SU(3) symmetry, and the 
(5(2j) = —1(5(2+) correlation to the quadrupole vibra¬ 
tional U(5) limit. Phenomenologically, the anharmonic 
vibration may also provide the <5(2)'’) = — (5(2^) corre¬ 
lation, although its spectral behavior is not observed in 
the shell-model TBRE. 

In particular, the invariance of under orthogonal trans¬ 
formation intrinsically provides the shell-model TBRE 
more opportunity to restore part of SU(3) properties, i.e. 
Q{2i) = ±(5(2J) correlations, even though these Q cor¬ 
relations are insensitive to the SU(3) rotational spectrum 
as expected based on the experimental survey [23|. On 
the other hand, IBMl Q{2^) = ±(5(2^) correlations al¬ 
ways favor low-lying rotational spectra, which indicates 
that the IBM is more strongly governed by the dynamic 
symmetry. The SU(3) group reduction rule also qual¬ 
itatively explains why the Shell Model more obviously 
favors the Qi^f) = —(5(2j) correlation compared with 
the IBMl. 

Low-lying Q correlations represent intrinsic nuclear 
collectivity, and thus are more sensitive to the wave- 
function detail than the spectrum. Therefore, the nuclear 
quadrupole collectivity may maintain in a far more deep 
level than the common realization based on the orderly 
spectral pattern. 
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